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Abstract 

We calculate secondary pairing in a model of a color superconductor with a quadratic 
gapless dispersion relation for the quasiquarks of the primary pairing. Our model mim- 
ics the physics of the sector of blue up and red strange quarks in gapless color-flavor- 
locked quark matter. The secondary pairing opens up a gap A s in the quark spec- 
trum, and we confirm Hong's prediction that in typical secondary channels A s oc G 2 S 
for coupling strength G s . This shows that the large density of states of the quadrati- 
cally gapless mode greatly enhances the secondary pairing over the standard BCS result 
A oc exp(— const /G). In all of the secondary channels that we analyzed we find that the 
secondary gap, even with this enhancement, is from ten to hundreds of times smaller 
than the primary gap at reasonable values of the secondary coupling, indicating that 
secondary pairing does not generically resolve the magnetic instability of the gapless 
phase. 



1 Introduction 



The exploration of the phase diagram of matter at ultra-high temperature or density is an 
area of great interest and activity, both on the experimental and theoretical fronts. Heavy- 
ion colliders such as the SPS at CERN and RHIC at Brookhaven have probed the high- 
temperature region, searching for the transition to deconfined quark matter. In this paper 
we discuss a puzzle that arises in a different part of the phase diagram, at low temperature 
but ultra-high density. Here there are as yet no experimental constraints, but calculations 
show that at sufficiently high density, the favored phase is color-flavor-locked (CFL) color- 
superconducting quark matter pQ (for reviews, see Ref. 0). 

The puzzle concerns the identity of the next phase down in density. Recent work [3! 
suggests that when the density drops low enough so that the mass of the strange quark can 
no longer be neglected, there is a continuous phase transition from the CFL phase to a new 
gapless CFL (gCFL) phase, which could lead to observable consequences if it occurred in the 
cores of neutron stars [3] . However, it now appears that some of the gluons in the gCFL phase 
have imaginary Meissner masses, indicating an instability towards an unknown lower-energy 
phase 0E1I3IE1. 

The gCFL phase is named after its most striking characteristic: the presence of gapless 
modes in the spectrum of quark excitations above the color-superconducting ground state. 
These include a gapless mode with an approximately quadratic dispersion relation E(p) oc 
(p — Pf) 2 , as well as gapless modes with the more typical linear dispersion relation E(p) oc 
\p — pf\. It seems likely that the presence of these modes is related to the instability 1 . It 
has therefore been suggested that one way to resolve the instability would be for the gapless 
quasiparticles to pair with each other ("secondary pairing"), leaving no gapless modes at all. 
In particular, it has been argued [TU] that because the quadratic gapless mode has so much 
phase space at low energy (its density of states diverges as E' 1 ^ 2 ) the secondary pairing will 
be greatly enhanced over the standard BCS value, which is based on pairing of modes that 
are linearly gapless. This offers the prospect that channels whose attraction seemed negligibly 
weak could become important once the primary pairing has created a quadratically gapless 
quasiquark, and these channels could then supply the required secondary pairing. It must 
be remembered, however, that in order to fully resolve the instability the secondary pairing 
gap parameter A s must be comparable to the primary pairing A p . If A s were significantly 
smaller then there would be a temperature range A s C T < A p in which there was primary 
pairing but no secondary pairing, and at those temperatures the instability problem would 
arise again. 

In this paper we investigate secondary pairing in a simple two-species model with a Nambu- 
Jona-Lasinio (NJL) interaction. We previously used a similar model to study the effects 
of gapless modes on photon and gluon screening masses jU]. Our model is essentially just 
one sector of the gCFL pairing pattern (the blue- up/red-strange sector), which is where a 
quadratically gapless quasiquark emerges, after primary pairing between the two species, in 
the Dirac C75 channel. 

The Cooper pairs must have an overall antisymmetric wavefunction, and because the 
secondary pairing is symmetric in color and flavor (pairing a given species with itself), its 

Note, however, that the instability can occur in fully gapped superconductors [S], for example at finite 
temperature UJ. 
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Dirac structure must be antisymmetric. We also restrict ourselves to rotationally invariant 
pairing, which leaves three possible Dirac structures for the secondary pairing: C, C75, and 
C7075. We give a detailed treatment of secondary pairing in the C7075 channel, because it is 
the only one that, for a single color and single flavor, is predicted to be attractive under an 
NJL interaction based on single gluon exchange (Ref. last 4 lines of Table 2). We work 
at zero temperature throughout. 

Our results confirm the conclusion of Ref. ^U] concerning the parametric form of the 
enhancement of the secondary pairing when it operates on quadratically gapless quasiparticles. 
We find that for a coupling strength G s in the secondary channel, A s oc G 2 S as compared to 
the BCS result A s oc exp(— C/G s ). This was predicted by Hong [TU] (note that our G s should 
be identified with Hong's effective interaction strength k, not with his This result can 

also be understood by a simple argument in the NJL model ^2] (see end of section |HJ). 

2 Two-species pairing formalism 

Our model contains two species of quark with primary pairing leading to anomalous self-energy 

{l^aClMlPI = \(<?l)ab ■ (2.1) 

The secondary pairing that we will discuss in most detail is the Dirac C7075 channel, 

i^aClQlMlPI = AA& • (2.2) 

The two-dimensional species space, indexed by a, b = 1,2, corresponds to the blue-up/red- 
strange sector of the combined color-flavor space in full QCD. As in the gCFL phase, the 
primary pairing is between species 1 and species 2, with the form (1,2) + (2, 1) given by the 
Pauli matrix o\. It is symmetric in the color-flavor indices (because it is antisymmetric in 
both color and flavor) and antisymmetric in the Dirac indices. It has spin zero (no spatial 
indices). The secondary pairing pairs each species with itself, and is also antisymmetric in 
the Dirac indices, and has spin zero. However it vanishes in the limit of massless quarks jU] 
because it pairs a left-handed quark with a right-handed quark. (For massless quarks with 
equal and opposite momenta this corresponds to pairing quarks with parallel spins, which 
would have to give a spin-1 state.) This means that C7075 secondary pairing can only occur 
if at least one of the flavors is massive. Fortunately that is the situation we are interested in, 
since one of our quarks is light (blue up) and the other is heavier (red strange). 

One might ask whether we have missed any important physics by using a model that only 
represents one sector of the gCFL pairing pattern. For example, there might be secondary 
pairing in other sectors that could feed back into the gap equations in our sector. However, as 
we will see below, even within our sector the feedback of the secondary pairing on the primary 
pairing is negligible, because the secondary pairing is so small, so there is no reason to expect 
significant contamination from other sectors. 

To treat the quark-quark condensation, which violates fermion number and allows quarks 
to turn into antiquarks, we use Nambu-Gor'kov spinors, which incorporate particles and an- 
tiparticles into the same spinor, (ip,ip T ). Our fermion fields \ are therefore 16-dimensional, 
arising from a tensor product of the 4-dimensional Dirac space, the 2-dimensional color-flavor 
space, and the Nambu-Gor'kov doubling. 
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The action for the quarks is 



A = \J x(pVS~\p)x(p)^. (2.3) 

The inverse propagator S 1 " 1 ^) is a 16 by 16 matrix: 
S -1 (po,p) = 

{PvloY + ft) ® 1 - 5fi ® ct 3 - 70 ® M ApC7 5 ® ai + A s C7o7s ® 1 
(A P C75 ® cti + A s C7o7 5 ® 1) T [(Pi/ToT"-/*) ® 1 + ® ^3+7o ® M] T 



(2.4) 



where in each tensor product the first factor lives in the Dirac space, and the second factor in 
the color-flavor space. The two species have an average chemical potential ft, and a chemical 
potential splitting 5fi which corresponds to the color and electrostatic potentials that enforce 
neutrality in the gCFL phase. The quark mass matrix is 

" = (2 ' 5) 

We can obtain the 4 branches (each 4-fold degenerate) of the full quark dispersion relation, 
including the effects of both primary and secondary pairing, by finding the values of the energy 
e(p) at which there is a pole in the full propagator 

det5 _1 (e(p),p) = . (2.6) 

If we index these solutions by a label a = 1 ... 16 then the grand canonical potential (or, 
loosely, the free energy) is 



f A I 

Q,(m u ,m s ,ft,6ft,A p ,A s ) = - / ^ e tt (p, m u , m s , Jl, Sfj,, A p , A s 



d 3 p Al A 
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(2tt)3 ^ G p + G. • (2 - 7) 



The effective couplings G p and G s are determined by details of the NJL model interaction. 
In this paper we will simply treat them as parameters. Our aim is to see how the secondary 
pairing A s depends on the coupling G s in the secondary channel. We do not attempt to use 
"realistic" values for G s (as might be obtained from commonly used NJL models of QCD) 
because we will be able to show that for any reasonable value of the secondary coupling (i.e 
f° r Gs \Gp) the secondary gap is too small to generically resolve the magnetic instability 
of a gap less phase. 

Our procedure is as follows. 

1. Choose values of ft, m u , and m s that are appropriate for quark matter in the core of 
a compact star, and an arbitrary ultra-violet cutoff A, significantly larger than ft. We 
used ft = 500 MeV with m u = to 100 MeV, m s = 160 to 300 MeV, and cutoff A = 800 
or 1000 MeV. 

2. Choose the desired gap parameter A p for the primary pairing. We varied A p from 
25 MeV to 75 MeV. To obtain a given value of A p one must choose a value of the 
coupling G p so that the free energy has a minimum at that value of A p . As we vary G p 
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we must also vary Sfx so that quasiparticle dispersion relation after primary pairing is 
always quadratically gapless. To simplify this process we first do it at m u = m s = 0, 
where we simply have to set 8/jL = A p to obtain quadratically gapless dispersion relations. 
I.e., we solve the gap equation 



2A p d 



G p dA p 



J 22\ e <x(P> m u = Q,m s = 0,/i, £/i = A p , A p , A s = 0) -^—^ (2.8) 



to obtain G p (m s = 0). We then turn on the desired value of m s , and retune 5fi to 
obtain quadratically gapless dispersion relations. It turns out that if we now use the 
primary gap equation (including the non-zero m s ) to determine what value of G p gives 
the desired primary gap parameter A p , the result is only slightly different (by < 10%, 
typically) from the value of G p {m s = 0). So in our results we actually use G p {m s = 0). 
Note that any "error" in G p is really just a rescaling of f3 in the secondary pairing gap 
equation (J2.2J) . whose only effect would be to slightly shift the line in Fig. ^ 

3. Study how the secondary pairing depends on G s . Recall that G p and G s both come 
from the same underlying NJL interaction with some coupling G. They arise from Fierz 
rearrangements of that interaction, so they are both of order G, with different numerical 
coefficients. It is therefore natural to define the ratio of the secondary to the primary 
effective coupling strength 

P = GJG P , (2.9) 

and to study how A s depends on f3. We expect /3 < 1 since the secondary pairing is 
by definition weaker than the primary pairing. For the C7075 secondary channel in a 
single-gluon-based NJL interaction, (3 ~ \ Our final step is therefore to solve the 
gap equation for the secondary pairing strength A s as a function of j3, 

d 3 p 



2A S d 
W P ~dA< 



/ I d p 

22\e a (p,m u ,m s ,p,,5fj,, A p , A s ) -^-^ (2.10) 



In principle one should solve coupled gap equations (J2.8J) and (J2.10J) . but we assume that 
A s <^ A p so, as we see explicitly below, the back-reaction of the secondary pairing on 
the primary gap is negligible. 



3 Results for the C7075 channel 

We present detailed results for the case fi = 500 MeV, m u = 0, m s = 300 MeV, with cutoff 
A = 800 MeV. For a primary coupling G = 7.625 x 10~ 6 MeV" 2 , which gives A p = 75 MeV, 
we tuned 5fi to 5/i qua( j = —26.54544010627093 MeV to obtain quasiquark dispersion relations 
that were quadratically gapless to within 10~ 13 MeV. We then calculated the secondary 
pairing as a function of (3. To confirm that the secondary pairing has negligible effect on the 
primary gap equation, we took the strongest secondary channel coupling, (3 = 1, which gave 
A s = 0.555 MeV, and re-solved the primary gap equation, including this value of A s in the 
primary gap equation. This led to a shift 5A P = 0.003 MeV, which is certainly negligible 
relative to A p = 75 MeV. As well as the gapless case, we also studied a slightly less negative 
value of 5fi, for which the quasiparticles had a small gap, and a slightly more negative value, 
for which there were two linearly gapless points. The results are plotted in Fig. ^ 
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Cy y 5 secondary pairing vs coupling 



A s /A p 




P=G S /G p 



Figure 1: The ratio of the secondary to primary pairing gap parameters A s /A p , as a function of 
the ratio of the coupling strengths in the primary and secondary channels (3 = G s /G p . In realistic 
QCD-like interactions we expect j3 < ^> an d our original assumption that A s -C A p turns out to 
be valid for all (3 in this range. For primary pairing with a quadratically gapless dispersion relation 
(squares) our results fit A s oc (3 2 , as predicted in Ref. JO] • For primary pairing with linearly gapless 
dispersion relation (small diamonds) we recover the BCS relation A s oc exp(— C / (3) at low [3. If the 
primary pairing leaves a gap (circles) then secondary pairing disappears at low coupling. 

1. Primary pairing giving exactly quadratically gapless quasiquarks (squares in Fig. [TJ) . 

This occurs when 5ji = 5/i qua d- Our NJL calculations of the secondary pairing follow 
the expected result, A s / A p = A(3 2 (our fit, straight solid line, is A = 5.287 x 10~ 3 ), all 
the way down to the lowest secondary couplings that we could probe. Note, however, 
that in the physically relevant range of couplings, (3 < 1, the secondary pairing is still 
suppressed relative to the primary pairing by a factor of 100 or more. (On the straight 
line fit, A a /A p (/3 = 1) = A = 5.287 x 10~ 3 ). 

2. Primary pairing giving gapped quasiquarks (circles in Fig. Q). 

We plot the case 5fi = —26.5454300, for which the gap in the spectrum is E gav w 
10~ 5 MeV. In this case the quasiquark spectrum looks quadratic at higher energies than 
this (see dashed line in schematic plot of dispersion relations, Fig. |2J) and so we expect 
the secondary pairing will only "notice" the gap if A s < 10~ 5 MeV i.e. A s /A p < 10~ 7 , 
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which corresponds to secondary coupling (3 ~ 0.01. The results of the explicit NJL 
calculation confirm this. For (3 > 0.01 we obtain a quadratic dependence A s /A p oc (3 2 , 
since the primary quasiquarks look quadratically gapless. At (3 ~ 0.01 the secondary 
pairing falls to zero, since there are no quasiquark modes within A s of the Fermi surface. 

3. Primary pairing giving linearly gapless quasiquarks (diamonds in Fig. QJ. 

We plot the case S/j, = —26.5454402, for which the height of the "bounce" in the pri- 
mary quasiquark dispersion relation is -Enounce ~ 10~ 7 MeV. The quasiquark spectrum 
looks quadratic at higher energies than this (see dash-dotted line in schematic plot of 
dispersion relations, Fig. |2J), and so we expect the secondary pairing will only "notice" 
the linearity if A s < 10 -7 MeV, i.e. A s /A p < 10 -9 , which corresponds to secondary 
coupling (3 ~ 0.0003. The results of the explicit NJL calculation confirm this. For 
(3 > 0.0003 we obtain A s /A p oc /3 2 , since the primary quasiquarks look quadratically 
gapless at energies greater that A s . At (3 ~ 0.0003 the secondary pairing changes to the 
BCS form, A s /A p = Aexp{-C/f3) (our fit is A = 3.333 x 10~ 8 , C = 0.001646), since 
the quasiquark modes within A s of the Fermi surface have linear dispersion relations, 
like those of unpaired fermions. 

We have checked that the results described above are generic. We have repeated the 
calculation for values of m s from 160 to 300 MeV, m u from to 100 MeV, and A p from 25 to 
75 MeV, for two different values of the cutoff A = 800 and 1000 MeV. In every case we found 
A s /A p < 0.01 at (3 = 1. 
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Figure 2: Schematic illustration of the dispersion relations of the quasiquarks (before secondary 
pairing) that underly the three sets of data plotted in Fig. ^ At higher energy scales, all three 
dispersion relations look quadratically gapless (upper panel). At energies < E g3iP it is possible to 
distinguish that one case is gapped (dashed line in lower panel). At energies < -^bounce it is possible 
to distinguish that the other has two linearly gapless points (dash-dotted line in lower panel). 
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There is a a simple calculation in the NJL formalism, pointed out by K. Rajagopal [T2] . 
which shows the origin of Hong's predicted A s oc (3 2 behavior for quadratically gapless quasi- 
quarks, and moreover enables us to see exactly how the C7075 gap is suppressed as m s — > 0. 
We first write down the generic form of the gap equation, 



1 = G : P 2 d P . (3.1) 

J y/E(p)* + A 2 

In the typical case where the quarks before pairing have linearly gapless dispersion relations 
E(p) oc \p — fj,\, the integral diverges logarithmically as A — > 0, and we obtain the BCS result 

1 « G / u 2 log(A/A) A « Aexp(-l/(GyU 2 )) . (3.2) 

In the case of a quadratically gapless dispersion relation, E(p) oc (p — fi) 2 / A p , if we assume 
(without justification, at this point) that the secondary pairing gap equation also takes the 
form (|3.1|) . with A replaced by A s , then the integral will diverge as an inverse square root of 
A s , so the gap equation for secondary pairing would become (using (J23 



1 oc G s fi 2 ^/A p /A s A s /A p oc ^G 2 p (3 2 . (3.3) 

which shows the A s oc (3 2 behavior. 

The full expression for the integrand on the right-hand-side of the gap equation in the C7075 
channel is very complicated, and we have not been able to find analytic approximations that 
would reduce it to a simple form that could be compared with (|3.1|) . However, the integrand is 
dominated by the contribution from the lowest-energy quasiparticle in the range of momenta 
where it would become quadratically gapless at A s = 0, and we have numerically evaluated 
that contribution in the case m u = 0, and found that it is well approximated (to within about 
10% for the values of //, m s , A p studied in this paper) by 

mU(p-p F ) 4 . (3 4) 



111 2 V Ap A* 2 

If we use this expression in the gap equation, we find 

1/2 / \ 3/2 



=^J i^j) => A s /A p ocm^G 2 p (3 2 



(3.5) 



This shows both the (3 2 dependence, and the suppression of the C7075 secondary pairing 
in the limit m s — > 0. We have checked the m s dependence by varying m s in our numerical 
calculations, and found that it is well-described by Eq. (J3.5j) . 

For secondary pairing in the other channels that are not suppressed in the m s — > limit 
(such as C75 [TH]) the factors of m 2 //i 2 in ()3.4|) would be absent and there would be a prefactor 
of /i 4 instead of rn^fi in (|3.5|) . For C7075 self-pairing of a single isolated flavor, as considered in 
Ref. [TT]). where the unpaired quasiparticles have linear dispersion relations, the gap integral 
has a BCS divergence. In that situation the factors of m 2 /fi 2 are exponentiated to give a 
much more severe suppression of the gap as m s — > 0, as seen in Fig. 4 of Ref. [TT] . 
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4 Other channels 



The detailed results presented above are all for condensation of the form (ipaC^Q^ipbjiPi = 
A s S a b, where the single-species pairing has the same sign for both species, (1, 1) + (2, 2). One 
could also construct secondary pairing of the form (i^aCjolb^b) xpi = A s (<r 3 ) a b where the two 
species each self-pair with opposite sign, (1, 1) — (2, 2). In Ref. [TT] these possibilities were not 
distinguished because the different colors and flavors were independent of each other. In the 
current context, however, there is primary pairing in the background, which connects the two 
species, so that the quasiquarks that undergo secondary pairing are not species eigenstates, 
but a momentum-dependent superposition of 1 and 2. This means that the relative sign of the 
secondary (1, 1) and (2,2) pairing is physically important. We have studied the C7075 (cr 3 ) a6 
channel, and we find essentially the same behavior as in the C7075 S a b channel: A s /A p < 0.01 
for < \ (first line of tabled). 

We have also analyzed the other Dirac-antisymmetric channels that are rotational scalars, 
arbitrarily assigning positive (attractive) interaction strengths G s to these channels. The 
results are summarized in table ^ For each channel we chose the same parameter values as 
used in Fig. □ namely p, = 500 MeV, m s = 300 MeV, A p = 75 MeV, m u = 0, A = 800 MeV, 
and we tuned 5fi to obtain a quadratically gapless quasiquark. In most channels we obtained 
results similar to those presented in Fig. ^ as the interaction strength drops from 1, A s /A p 
drops as (3 2 . The only exception is the CS a b channel, where we find that there is no secondary 
pairing for ^ 1. However, among the channels studied we find significant variation in the 
strength of the secondary pairing. The potentially physically relevant region is where the 
secondary channel interaction strength G s is smaller than the primary channel interaction 
strength G p . In table ^ we therefore show A s /A p at (3 = ^ and 1. We see that in the 
C^f^{a 3 ) a b and C(a 3 ) a b channels the secondary pairing is about ten times weaker than the 
primary pairing when (3 = |, and that in order for A s to be of the same order as A p , the 
secondary channel would have to be as strongly attractive as the primary channel. 



Channel 


A s /A p 


Channel 


A s /A p 




(3=\ (3 = 1 




(3 = \ (3=1 




0.0015 0.0074 


C7o75 (c^U 


0.0024 0.012 


C75 $ab 


0.00019 0.026 


C75 (cr 3 ) ab 


0.093 0.72 


C Sab 


< 10- 11 < 10- 11 


C (0-3)06 


0.067 0.53 



Table 1: Secondary pairing in various rotationally- invariant single-species channels. Param- 
eters were the same as in section 01 The secondary channel would have to be as strongly 
attractive as the primary channel ((3 = G s /G p = 1) for the secondary pairing to be compa- 
rable to the primary pairing. As (3 drops below 1, A s rapidly drops to a small fraction of 
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5 Conclusion 



Our results validate the intuition that when the quark dispersion relations (including primary 
pairing) are quadratically gapless, the very large density of states near the Fermi surface 
should lead to a parametric enhancement of the secondary pairing, relative to BCS pairing, in 
weakly coupled channels. However, it is clear that even with this enhancement, the secondary 
pairing is not typically of comparable magnitude to the primary pairing. 

For the case of the C7075 channel that is predicted to be attractive in single-gluon- 
exchange-based NJL models, the secondary pairing is at least 100 times weaker than the 
primary pairing. The other rotationally-invariant secondary channels are predicted to be re- 
pulsive (G s < 0), but we have studied how they would behave if they were attractive. We 
find that for G s /G p ^ |, they are at least ten times weaker than the primary pairing. For 
G s /Gp approaching 1 the secondary pairing could become important in certain channels, but 
that would essentially correspond to assuming that they were not secondary after all, and 
the competition between the "secondary" and primary channels would have to be taken into 
account. We therefore conclude that in spite of the great enhancement provided by the large 
density of states of a quadratically gapless quasiquark, secondary pairing will not generically 
resolve the magnetic instability of the gapless phase. 

This result was obtained under the most favorable conditions for secondary pairing (quadratic 
gapless dispersion relation), so we expect that it is robust. For example, in full gCFL matter 
it might turn out that secondary pairing affects the neutrality condition that it is no longer 
physically correct to tune Sfi in the blue-strange/red- up sector so as to obtain a quadratic dis- 
persion relation. However, as is clear from Fig. Q this can only further suppress the secondary 
pairing gap. This means that our conclusion will also apply to the (blue-down/green-strange) 
sector of the gCFL pairing pattern, whose quasiquark is linearly gapless with an energy scale 
-^bounce ~ A*e/2 which is typically tens of MeV. Our results for linearly gapless quasiquarks 
(diamonds in Fig. show that weakly attractive secondary channels could readily open up 
a tiny (BCS-suppressed) gap at the linearly gapless points. However, in order to resolve 
the Meissner instability we would need secondary pairing of the same order as the primary 
pairing in this sector as well as in the blue-up/red-strange sector. When A s > -Enounce the 
quasiquarks will behave as if they were quadratically gapless, so our analysis also applies to 
secondary pairing in the blue-down/green-strange sector, and our conclusions apply there as 
well. 
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